Resonant Shattering of Neutron Star Crusts 
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The resonant excitation of neutron star (NS) modes by tides is investigated as a source of short 
gamma-ray burst (sGRB) precursors. We find that the driving of a crust-core interface mode can 
lead to shattering of the NS crust, liberating ~ 10''^ — 10*^ erg of energy seconds before the merger 
of a NS-NS or NS-black hole binary. Such properties are consistent with Swift/BAT detections of 
sGRB precursors, and we use the timing of the observed precursors to place weak constraints on 
the crust equation of state. We describe how a larger sample of precursor detections could be used 
alongside coincident gravitational wave detections of the inspiral by Advanced LIGO class detectors 
to probe the NS structure. These two types of observations nicely complement one another, since 
the former constrains the equation of state and structure near the crust-core boundary, while the 
latter is more sensitive to the core equation of state. 

PACS numbers: 97.60.Jd, 97.80.-d, 98.70.Rz, 95.85.Sz 



Introduction. Short-hard gamma-ray bursts (sGRBs) 
are intense bursts of gamma-rays with isotropic equiv- 
alent energies of E^so ~ 10^" — 10^^ erg, and durations 
< 2 s [T] . The leading model for the progenitors of sGRBs 
is the merger of a neutron star (NS) with another NS or 
the tidal disruption of a NS by a black hole (BH) [2 [3]. 
Studying the response of NSs in dynamic systems pro- 
vides a window into the physics of their crust and core, 
probing material up to many times nuclear density in 
states inaccessible to terrestrial laboratories. 

Flares occurring before the main GRB — precursor 
flares — were recently identified for a handful of sGRBs 
[H [5]. In particular, precursors ~ 1-10 s prior to the 
main flare were detected with high significance for three 
sGRBs out of the 49 considered in Troja et al. [S]. Until 
these discoveries, precursor flares had only been identi- 
fied prior to some Long Gamma Ray Bursts (e.g. [SIT]). 
Detailed sGRB precursor mechanisms remain largely un- 
explored, but since the precursors occur when the NSs 
are strongly interacting, they present an exciting oppor- 
tunity to learn about the interior or region around the 
NSs. Magnetospheric interaction has been proposed as 
one explanation for the precursors, but is unlikely to be 
sufficiently strong unless the magnetic fields exceed mag- 
netar strength [S]. Direct tidal crust cracking has also 
been suggested [21 [TO] • The crust cracks when the crust 
breaking strain et — 0.1 [TT] is exceeded by the direct 
tidal ellipsoidal deformation [T^] 
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where hi is the quadrupole shape Love number of the 
deformed body, Mia = M/IA Mq and R12 = R/12 
km are its mass and radius, respectively, q is the mass 



ratio of the binary, and fgw is the gravitational wave 
frequency (twice the orbital frequency); post-Newtonian 
corrections to the tidal potential (PN) are found in 13J. 
At the frequency when SR/R ^ eb, the gravitational wave 
inspiral timescale [14j . is 
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where M = M^^^M^^^/{Mi + JVh)^/^ is the chirp mass. 
Since this timescale is so short, unless the GRB emission 
is delayed for seconds after the objects coalesce (not gen- 
erally expected), direct tidal deformation cannot explain 
the majority of observed precursors. 

In this letter we present an alternate tidal mechanism 
that can shatter the crust and cause fiares: the excita- 
tion of a resonant mode by periodic tidal deformation. 
We identify a particular mode that has an amplitude 
concentrated at the crust-core interface, and show that 
its coupling to the tidal field is sufficient to fracture the 
crust. Further driving of this mode can continue to cause 
fractures and deposit energy into seismic oscillations un- 
til the elastic limit is reached, releasing the entire elastic 
energy of the crust as it shatters. Such an event may be 
observable as a precursor flare. 

Steiner and Watts [TS] showed that the identification of 
quasiperiodic oscillations in magnetar flares with toroidal 
shear modes can constrain physical parameters such as 
nuclear symmetry energy and we similarly demonstrate 
how multi-messenger probes of sGRB precursor bursts 
can be used to measure the interface mode frequency and 
constrain the crust equations of state. 

Analysis. The tidal excitation of NS normal modes 
in a binary has been well studied for fluid modes (e.g. 
[TOl [TT]). and Reisenegger and Goldreich [TBI hinted at 
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possible associated electromagnetic signatures. Including 
the solid crust of a NS into the mode analysis changes the 
normal mode spectrum significantly [19]. In addition to 
the fluid modes of the core (e.g. f-modes, g- modes), the 
solid crust induces several new modes, including crustal 
shear modes and interface modes between the crust and 
core. 

For the present study we focus attention on the mode 
most likely to be tidally excited and subsequently crack 
the crust. The toroidal crust shear modes have large 
shear strains, but do not couple significantly to the 
spheroidal tidal field and are negligibly excited. Low or- 
der core g-modes can couple to the tidal field [16], but 
they do not penetrate the crust. The f-mode of the NS 
has the strongest tidal coupling, and can strongly perturb 
the crust. However, the f-mode frequencies are typically 
J> 10'^ Hz. These frequencies are not reached by the bi- 
nary before merger, or are in resonance only at late times 
(see eq. [2]). This would make any observable signal of 
this mode difficult to separate from the sGRB emission. 
In contrast to these other modes, the I = 2 spheroidal 
crust-core interface mode (first identified by |19)V or i- 
mode, is ideal for resonant crust fracture as it couples 
significantly to the tidal field and has a large shear strain 
near the base of the crust. In addition, its frequency is 
low enough (~ 100 Hz) for tidal resonance to be well sep- 
arated from the merger (tg^ ~ 2 s at fgw — 100 Hz using 
eq. 1). 

We construct background models of spherical neutron 
stars with the Oppenheimer-Volkoff equations. For the 
discussion below, we use the IAMq (12km) NS with 
Skyrme Lyon (SLy4) EOS [I^ and a core/crust tran- 
sition at baryon density of rib ~ 0.065 fm~^ as a fidu- 
cial case. Taking the NS perturbation equations from 
[l8[ [To] we solve for the crust-core i-mode with a crust 
shear modulus given by (2U] (which ignores free neutron 
interactions, and superfluid entrainment effects), 

^ 0.1194 n,{Zef 

l + 0.595(173/r)2 a ' ^ ' 

where Z is the atomic number of the ions, rii is the ion 
density, F = {ZeY /akbT is the ratio of the Coulomb 
to thermal energies and a = (3/47mi)^/^ is the average 
inter-ion spacing. We assume T w 10* K and ignore any 
magnetic field since a field strength ^ lO^^G is required 
for the magnetic energy density to be comparable to the 
shear modulus fx ~ 10'^° ergcm"^ at the base of the crust, 
and such a field has decay time that is likely much shorter 
than the binary lifetime of the system [3T]. The exact 
mode frequency and eigenfunction depend on properties 
near the crust-core boundary including EOS, crust shear 
modulus and crust-core transition density. There is also 
a weaker dependence on bulk properties of the star such 
as its mass. For this model the i-mode frequency is found 
to be /mode — 188 Hz, with eigenfunction shown in Fig. 
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FIG. 1: Left Panels: The I = 2 crust-core i-mode eigenfunc- 
tions as a function of radius where the Lagrangian displace- 
ment is i{r,e,<i}) = U{r)Y2±2{e,<t))f + V {r)VY 2±2{e , <f)) ■ A 
1.4M0 12 km NS with SLy4 [15| equation of state is used for 
the background model. Right Panel: Cross section of the NS 
deformation due the mode. The binary companion is located 
along the horizontal axis, about which the deformation has 
rotational symmetry. The fluid core is ellipsoidally deformed 
towards the binary companion. The solid crust is sheared 
strongly and compressed to make up for the deformation of 
the core such that the outer surface is only slightly perturbed 
radially. The horizontal strain and radial deformation peak 
at the base of the crust, making excitation of this mode ideal 
to crack the crust. The solid black line denotes the surface of 
the unperturbed star, while the dashed black line denotes the 
unperturbed crust/core transition radius. 

As the binary inspirals, it sweeps through successively 
higher frequencies until it passes through a resonance 
with the mode frequency. The timescale of during which 
the tidal driving is approximately phase coherent with 
the i-mode is estimated with a random phase approxi- 
mation to find tics ^ ^/tgw/T^fgw [K] , 

where we have used the fact that /mode ~ fgw during 
the resonance. During this time, we estimate the degree 
of resonant excitation of the i-mode using the formalism 
of Lai [TS]. From this we calculate the dimensionless 
overlap integral between the I = 2 i-mode and the tidal 
field finding 

Q = J d\p e ■ V[r2Y2±2(f?, 0)] ^ 0.041 , (5) 

where | = U{r,t)Y2±2'r + V{r,t)'VY2±2 is Lagrangian 
displacement mode eigenvector normalized such that 
/ d?xp$, ■ i* = MR"^, and Y2±2{S,(t>) is the spherical 
harmonic with I = 2 and m = ±2 with both m modes ex- 
cited equally, and the numerical factor on the right-hand 
side is from our calculation using the SLy4 equation of 
state. 

During the resonance, the energy of the i-mode in- 
creases quickly over the timescale t^cs- The maximum 
possible energy the mode can attain was explored by |16) , 
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ignoring crust fracture, orbital back reaction and other 
non-linear effects. Using equation (6.11) from this work, 
we estimate the maximum possible energy of the i-mode 
to be 

5/3 

i^^ax ^ 5 X 10^" erg fUsQloiM^.i" ^(2 9 ' ' ' 
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where /iss = /mode/188 Hz and Qom = Q/0.04. We 
next compare this to the mode energy needed to create 
a break in the crust E\^. For a given amplitude mode, we 
can numerically calculate the dimensionless strain e (see 
e.g. |22)). and identify the location of maximum strain 
within the crust. In general, this is found at the base 
of the crust and concentrated near the NS equator. We 
numerically find the minimum amplitude |^b| needed for 
e = eb at some location within the crust, which in turn 
implies a minimum energy. 
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where eo.i = fb/O.l. Since _Eb ^ E^-mx the i-mode will be 
amplified enough to break the crust early in the resonant 
energy transfer. 

These initial breaks have a characteristic size similar 
to the thickness of the crust Ar ~ O.li?, and thus release 
elastic energy ~ e^/xAr"^ ~ lO"*"^ erg from the excited 
mode. This energy is converted primarily into a broad 
spectrum of seismic waves, peaked at characteristic fre- 
quency ~ (^/yo)-^/^/(27rAr) ~ 200 Hz. Seismic waves at 
these low frequencies do not couple efficiently to the mag- 
netic field, so the energy remains mostly in the crust [23] . 
But the fracture is not limited to a single event; the crust 
continues to fracture as the location of maximum strain 
moves around the star, since the NS is not tidally locked 
[Mj and the crust heals extremely quickly ^llj. Seismic 
energy injected by these fractures can build up in the 
crust until it reaches the elastic limit and the crust shat- 
ters. The elastic limit can be estimated as 
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The crust fractures cause the i-mode amplitude to satu- 
rate at amplitude |^b|, with the tidal energy transfer rate 

m 
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where $3 = J2im the gravitational potential 

due to the companion. We have the components 
V$2±2 = (/gw/7r)2g(g + l)-^W2±2V{r^Y2±2), where 
W2±2 = (67r/5)i/2 and taking a^VSt = 27r/modci:, we 
can use equations ([5| and Q to rewrite this expression 
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TABLE I: Resonant mode properties for the / — 2 i-mode. 
The background star is taken to be a 1.4 Mq NS, with various 
equations of state given in [TS1. The crust/core transition 
baryon density is fixed to be nt = 0.065 fm~^ for each model. 



where E^^ = Eb/W^^ erg. This high rate im- 
plies that only a small fraction of the resonant time, 
-EoiasticZ-E'tidai ~ 10"^ s < tj-os, IS required for the en- 
ergy in the crust to build up to the elastic limit. 

When this limit is reached and the crust shatters, this 
scatters the elastic energy and mode energy to higher 
frequency oscillations, which can couple efficiently to 
the magnetic field j25|, I26j . If emission occurs over the 
observed precursor timescale this implies a luminosity 
L ^ {Eb + E,i^tic)/0-ls ~ 10*8 erg s'^. If the mag- 
netic field is strong enough that the perturbations are 
linear, B > 10^^ G RuiL/lO'^^ergs-^)^/'^ , the energy 
propagates as Alfven waves along open magnetic field 
lines, and the resultant emission will be non-thermal. For 
weaker fields the magnetic perturbations are highly non- 
linear and can generate strong electric fields, which can in 
turn accelerate particles to high energy, and spark a pair- 
photon fireball with a near-thermal spectrum (e.g. [3 ). 
The observed temperature of such a flare depends sensi- 
tively on the baryon load carried in the fireball ^27. ,2S\. 
The total fluence of the precursor fiare can also provide 
an estimate for a lower bound on the breaking strain, 
since both Ef, and £'oiastic scale with e^. 

Motivated by the above calculations, we next explore 
the range of energetics and timescales expected from a 
sample of NS crust EOSs, and describe what constraints 
the observations may provide between them. 

Model Exploration. Using the theoretical framework 
described above, we perform the same analysis with the 
EOSs APR, SkI6, SkO, Rs, and Gs, following Steiner and 
Watts ^15j and references therein. Results are summa- 
rized in Table in O ur main conclusion is that the energet- 
ics are robust, since they mainly depend on the density of 
the crust and displacement eigenfunction, which do not 
vary greatly between the models. Much like the toroidal 
shear modes [l5l|29], the i-mode frequency is only weakly 
dependent on NS mass, but varies with shear speed at the 
base of the crust. This implies that the resonant excita- 
tion of the i-mode would occur at different times during 
the merger event. 

To quantify the relationship between the crust EOS 
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and the time at which the i-mode resonance occurs, we 
plot the i-mode and gravitational wave frequency versus 
time until (PN) coalescence, tc — t — Stg^/S (see e.g. 
[3D]), in Figure [2j The dashed lines trace the leading 
order frequency evolution for a given chirp mass TW, go- 
ing from left to right. When the dashed line intersects 
a colored column, it indicates the time and frequency at 
which resonance occurs. From this set of EOSs and M, 
a wide range of timescale are possible, from < 0.1 s up to 

20 s before merger. Also plotted as horizontal dotted 
lines are the observed precursor times reported in [5]. 
Although this comparison assumes that the main flare 
is nearly coincident with the binary coalescence, certain 
constraints can already be inferred. The relatively high 
frequency of the i-mode for the SLy4 EOS means that 
the resonance only occurs at late times, close to merger. 
Only if < IMq can such a model give timescales sim- 
ilar to the shortest precursors and the longer precursors 
may be especially difficult for this model to replicate. 
Other EOS models, such as Gs, Rs, SkI6, and SkO, are 
largely consistent with the timescale of precursors, but as 
a larger sample of precursor observations are made, dia- 
grams such as this will be useful for constraining EOSs. 

A binary with unequal mass NSs may excite two pre- 
cursor flares separated by a small time delay, due to the 
slight difference in the i-mode frequency. However, the 
two precursors (13 s, 0.55 s) observed in GRB 090510, 
are too far separated to both be explained by our reso- 
nant shattering model of precursors, using two NSs with 
the same EOS. The 0.55 s flare may alternatively be ev- 
idence of direct crust cracking [TD] and a delayed main 
GRB burst, the formation of a hyper-massive magnetar 
before collapse into a black hole [31] , or some other flare 
mechanism. 

Discussion. We explored the resonant excitation by 
tides of a mode that is concentrated at the crust/core 
boundary of NSs. We demonstrated that the resonance 
occurs between ~ 0.1 — 20 s prior to merger in NS-NS 
or NH-NS binaries. Further work remains to be done 
exploring the details of this model, including the effects 
of damping on the mode excitation, the effect of more 
realistic NS structure, and the detailed physics of the 
magnetospheric emission. However, we have shown that 
the energetics of the release of mode and elastic energy 
and the timescale at which the resonance occurs are sug- 
gestive of the precursors of sGRBS. Using this theoretical 
framework we demonstrated that interesting constraints 
can be placed on the NS crust EOS with comparisons to 
precursor observations. 

The direct phase change of the gravitational wave- 
form due to the resonant excitation of the mode, A(j) ^ 
(tgwEh) / (toThitEorhit) ^ 10"'^ rad, is too small to be di- 
rectly measured for signal to noise (SNR) <^ 1000. How- 
ever, coincident timing between the 7-ray burst detectors 
and the GW detector would allow precise determination 
of the mode frequency, coalescence time, main burst de- 
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FIG. 2: The time until PN coalescence {tc — t) as a, function 
of gravitational wave frequency. The dashed lines show the 
frequency evolution of inspiraling binaries for different chirp 
masses M as labeled in units of Mq. A given binary moves 
from left to right in time. The colored columns show the 
resonance frequencies /mode = /gw of a set of crust EOSs 
from [15], over a neutron star mass range of 1.2 (higher 
frequency) to 1.7Mq (lower frequency). We take 1.2Mq as the 
smallest plausible companion mass, giving an upper bound on 
the precursor times for each EOS. NS-NS systems will have 
chirp masses of 1.0 — l.SM©, and NS-BH systems with 10 — 
2OM0 BH have chirp masses of 2.7 — 4.5M0. The precursor 
times for the GRBs reported in [5] are plotted as horizontal 
dotted lines. 



lay time, and chirp mass. With parameter extraction 
from the GW inspiral at the detection threshold with 
SNR ^ 10, the dominant error in determining the res- 
onant frequency is due to the uncertainty in the timing 
of the precursor flare, which is of order the precursor 
duration. This implies that the mode frequency can be 
determined to fractional accuracy A/// ^ O.ls/igw 
2% (A^/1.2)^/^/*gQ . Such a measurement would allow us 
to tightly constrain the NS physics and parameters that 
determine the mode frequency. This is complementary 
to the constraints given by GW coalescence measurement 
alone, which are sensitive primarily to the core EOS (e.g. 
[321133]). 

Resonant shattering precursor flares are likely to be 
fairly isotropic, and thus may be observable even for 
sGRBs where the main flare is beamed away from the 
Earth. Such flares may also be a source of electromag- 
netic emission for higher mass ratio, lower spin NS-BH 
mergers where the neutron star does not disrupt to pro- 
duce a torus and main sGRB flare [M] |3S] . 
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